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Abstract 

Our calculation of the total cross section for inclusive production of ti pairs in hadron collisions 
is summarized. The principal ingredient of this calculation is resummation of the universal 



. leading-logarithm effects of gluon radiation to all orders in the quantum chromodynamics 

CsJ ' couphng strength, restricted to the region of phase space that is manifestly perturbative. We 

^ ' present predictions of the physical cross section as a function of top quark mass in proton- 



antiproton reactions at center-of-mass energies of 1.8 and 2.0 TeV. 



1. Introduction 



In this report we summarize our calculation of the inclusive cross section for top quark 
production^. In hadron interactions at collider energies, the main production mechanisms for 
^ ! top-antitop quark (ti) pair production, as modeled in perturbative quantum chromodynamics 

• (pQCD), involve parton-parton collisions. The gluonic radiative corrections to the lowest-order 

channels create large enhancements of the partonic cross sections near the top pair production 
threshold!. The magnitude of the 0{al) corrections implies that fixed-order perturbation 
theory will not necessarily provide reliable quantitative predictions for ti pair production at 
Fermilab Tevatron energies. A resummation of the effects of gluon radiation to all orders in 
perturbation theory is called for in order to improve the reliability of the theory. US. 

We use the notation = a;s(yu)/7r, where /i is the common renormalization/factorization 
scale of the problem. Unless otherwise specified, a = a{fi = m) where m is the mass of the 
top quark. For the subprocess 

+ t{p,) + i{p2) + g{k), (1) 
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we use the partonic invariants! 

s = {ki + /c2)^, ti = {k2 - - m^, ui = (ki - - ■> (1 - ^)^'^ = s + ti+Ui. (2) 

Through next-to-leading order (NLO), keeping only the leading logarithmic contributions, we 
write the total partonic cross section in the MS scheme as 

ai°+'^(r/, m^) = dz{l + a2Cij ln'(l - z)]a[A7], z, m') , (3) 



z, m 



where r] = s/^Am'^) - 1, Cij = Cp = 4/3 {Qj = Ca = S) for ij = qq (gg), (T^(?7, ^, 
■j^o'fjiV: ^1 rn"^), and afj is the unpolarized Born partonic cross section. We invoke universality 
with the Drell-Yan case {II production). Because the finite-order leading logarithms are iden- 
tical in the ti and // cases, we can resum them in tt-production with the same function we use 
in the Drell-Yan case. We find 

aij{r],m^) = [ dzH{z, a)a-Ar], z,m'^) . (4) 

Jl-4(l+r;)+Vl+,7 



The kernel of the hard part isi 

/•ln(-^) 

niz, a) = l+ ' c/xe^(^'") Qj(^^ ") • (5) 

The kernel in Eq. depends solely on the resummation exponent E{x,a), either explicitly, 
or through the functions Qj which depend exclusively on E. 

2. The Resummation Exponent 

For the Drell-Yan process, the exponent in moment spaceH in the Principal Value Resum- 
mation (PVR) approach^ may be written in the MS factorization scheme as 

Fix, a) = -g^'^ / dC \ . / ^ ^T— ITTT • ^ 
Jp 1 — C l + ao2lnA 

Here P is a principal- value contour^ and g^^'' = 2Cij. It is important to note that Eq. (P) can 
include in general all large logarithmic structures in the Drell-Yan case. It has a perturbative 
asymptotic representationB 

N(t)+1 p+1 

Fix, a) ~ Fix, «, Nit)) = (7(1) E E Hp^' ' (7) 

p=l j=0 
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where 

sj,p = -&r'(-l)^+^2%+i_,(p - l)!/j!, (8) 

and r(l + z) = Y^'^^QCkz'' , where F is the Euler gamma function. This representation is vahd 
in the moment-space interval 

1 < X = lnn < t. (9) 

This range of vahdity has the consequence that terms in the exponent of the form In*^ n 
are of order unity, and terms with fewer powers of logarithms, a'^ In'^"'" n, are negligible. This 
explains why resummation is completed in a finite number of steps. In addition, we discard 
monomials a'^ In^ n in the exponent because of the restricted universality between the ti and 
U processes. The exponent we use in our calculations is the truncation 

Nit)+1 

E{x,a,N) = g^^^ ^ (10) 
p=i 

with the coefficients Sp = Sp+i,p = b2^^2P/p{p+ 1). The number of perturbative terms N{t) in 



Eq. is obtained by optimizing the asymptotic approximation 



E{x,a) - E{x,a,N{t)) 



minimum. An excellent numerical approximation is provided by the fit N{t) ~ [t — 3/2] , 
where the integer part is defined as the closest integer from either direction. Throughout this 
paper, we use the two-loop formula for the fixed coupling strength a{m). 

In Fig. la we illustrate the validity of the asymptotic approximation for a value of t 
corresponding to m = 175 GeV. Optimization works perfectly, with N{t) = 6, and the plot 
demonstrates the typical breakdown of the asymptotic approximation if were to increase 
beyond N{t). This is the exponential rise of the infrared (IR) renormalons, the (p — 1)! growth 
in the second term of Eq. (j^) . As long as n is in the interval of Eq. , all the members of the 
family in n are optimized at the same N{t), showing that the optimum number of perturbative 
terms is a function of t only. 

It is valuable to stress that we can derive the perturbative expressions, Eqs. (|^), (||) and (H), 
without the PVR prescription, although with less certitude. We begin with the unregularized 
form of Eq. (^), i.e., with the integral over x on the real axis. Expanding the inner integrand 
as a Taylor series around a, we find the same result as Eq. (|^, the only and major difference 
being that we do not know the asymptotic properties of this series in the full range of moments 
n. For a fixed t and n, one may use the monotonicity behavior of the corresponding partial 
sums to try to determine an upper limit for the number of terms in Eq. (|^. This procedure 
is illustrated in Fig. lb for m = 175 GeV. We note that beyond a certain range of A^, the 



3 




7 8 9 10 11 12 13 14 15 IS 17 18 

N, m-175 GeV 



2 3 4 5 6 7 



(a) 



(b) 



I 13 14 15 16 17 

N, M=175 GeV 



Fig. 1. Optimum number of perturbative terms in the exponent (a) with PVR (sohd family is 
for PVR, dashed for perturbative approximation, both famihes increasing, for parametric values 
n = 10, 20, 30, 40) and (b) using monotonicity of the partial sums (solid family for n e {10^, 10'^} in 
steps of 100, dotted for n = 2, 3 x 10^ and dashed for n = 1, 2, 3 x 10''). 

exponent increases factorially, a demonstration of both the asymptotic nature of the series 
and of the effect of the IR renormalons. A range of optimum N can be determined where 
the growth of the sum reaches a plateau, before the factorial growth sets in at large A^. The 
plateau is centered around Nopt — N{t) = [t — 3/2], in agreement with PVR. 

3. The Resummed cross section 

The general form of the PVR cross section is given by Eqs. (H) and @. The functions 
Qj{x, a) in Eq. (^) appear during the inversion of the Mellin transform^ n ^ z. Specification of 
the perturbative regime in momentum space follows from general expressions for the inversion 
of the Mellin transform and the meaning of the successive terms in this inversion, once their 
perturbative approximations are used. The functions Qo and Qi, can be calculatedB: 



Qo = - sin(7rPi)r(l + Pi), Qi ~ 21(1 + Pi)P2 cos(7rPi)*(l + Pi) 



where P^ = [d^ /k\d^x)E and ^ = \E''^°) and '^^^^ are the usual Polygamma functions. For 
simplicity we include in the expression for Qi only terms that generate corrections starting 
at 0{a). Qi contributes one less power of x than a in the integrand of Eq. (|^), and it is 
formally subleading relative to the contribution of Qq. Nevertheless, from Eq. ([TT|) we see 
that this suppression is not true for values of x such that Pi(a;, a) ^ 1. We conclude that the 
perturbative region in momentum space is defined by the inequality constraint 



Piixz^a) < 1 , = In 



l-z 



(12) 
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Since we intend to resum leading logarithms only, our main result for the perturbative re- 
summed partonic cross section, denoted by af^'^^^* , is 



cr, 



R;pert 



dz 



1+ r'cia;e^(^'")Pi(x,a) 



20 



l-4(l+r;)+Vl+»? 



(13) 



where zq is the end-point calculated from Eq. (|T2]). In order to achieve the best accuracy 
available we wish to include in our predictions as much as is known theoretically. Our "final" 
resummed partonic cross section can therefore be written^ 



pert/ 2 2\ 



0", 



R;pert 



I (0+1)/ 2 2 



(14) 



R;pert 



The second term is the part of the partonic cross section up to one-loop that is included in 
the resummation, while the last term is the exact one-loop cross section!. 

It is useful to translate our definition of the perturbative regime directly into a statement 
about the perturbative region in t]. Our perturbative resummation probes the threshold down 
to the point rj > r]o = {1 — zq)/2. Below this value, perturbation theory, resummed or 
otherwise, is not to be trusted. The physical cross section for each production channel is 
obtained through the factorization theorem. 



aij{S, 



S 
4m2 



Am 



S 



(15) 



where the parton flux is = ^fi/h^{x,jj?)fj/h2{ylx,iJ?). The total physical cross 

section is obtained after incoherent addition of the contributions from the the qq and gg 
production channels. We ignore the small contribution from the qg channel. We adopt 
one common scale /i for both the factorization and renormalization scales. A quantity of 

phenomenological interest is the differential cross section fHi(^£!ij!z) _ j^g integral over 7] is, 
of course, the total cross section. In Fig. 2 we plot these distributions for m = 175 GeV, 
\/S = 1.8 TeV and fi = m. We observe that, at the energy of the Tevatron, resummation is 
significant for the qq channel and less so for the gg channel. We show the total tt-production 
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Fig. 2. Differential cross section da/drj in the MS-scheme for the (a) qq and (b) gg channels: Born 
(dotted), NLO (dashed) and resummed (solid). 

cross section as a function of top mass in Fig. 3. The central value of our predictions 
is obtained with the choice /i/m = 1, and the lower and upper limits are the maximum and 
minimum of the cross section in the range of the hard scale fi/m E {0.5,2}. Our prediction 
of Fig. 3a is in agreement with the datai. We find a*\m = 175 GeV,v^ = 1.8 TeV) = 
5.52lo;42 pb. In Fig. 3b we present our predictions for an upgraded Tevatron operating 
at y/S = 2 TeV. Our cross section is larger than the NLO one by about 9%. We predict 
a"(m = 175 GeV, = 2 TeV) = T.^et^oil pb. Over the range fx/m E {0.5, 2}, the band of 
variation of the strong coupling strength as is a generous ±10% at m = 175 GeV. 



agg{m = 150 GeV; DIS) 


fi/m = 0.5 


fi/m = 1 


fi/m = 2 


NLO 


9.42 


9.31 


8.57 


pert 


9.76 


9.92 


9.31 


LSvN(/io = 0.1m) 


7.9 


10.0 


9.7 


agg{m = 150 GeV; MS) 


fi/m = 0.5 


fi/m = 1 


fi/m = 2 


NLO 


2.51 


2.22 


1.81 


peri 


2.53 


2.30 


1.89 


LSvN(/io = 0.2m) 


1.76 


4.38 





Table 1. Physical cross sections in pb at m = 150 GeV and ==1.8 TeV. The corresponding LSvN 
predictions are also shown, where /iq is their IR cutoff. 

Comparing with Laenen, Smith, and vanNeerven (LSvN)i, we find our central values are 
10 — 14% larger, and our estimated theoretical uncertainty is 9 — 10% compared with their 
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p t t+X, Vs=1.8 TeV 



p p^ t t+X, Vs=2.0 TeV 



140 150 160 170 laO 190 200 210 220 230 240 250 260 140 150 160 170 180 190 200 210 220 230 240 250 260 

m(GeV) m{GeV) 

Fig. 3. Inclusive cross section for top quark production in the MS-scheme. The dashed curves show 
our perturbative uncertainty band, while the solid curve is our central prediction: (a) VS = 1.8 TeV 
and (b) \/5 = 2 TeV. 



28% — 20%. Our Born cross section, however, is about 3 — 5% larger than the LSvN Born 
cross section due to the different parton distributions used in the two calculations, including 
differences in A which alone account for half or more of this increase. The two predictions have 
overlapping uncertainties and are, in this sense, in agreement. It is important to stress that 
the theoretical uncertainties are estimated in quite different ways in the two methods. We use 
the standard //-variation, whereas LSvN obtain their uncertainty primarily from variations of 
their undetermined IR cutoffs. One of the advantages of a resummation calculation should be 
diminished dependence of the cross section on /i, less variation than is present in fixed-order 
calculations. Table 1 shows that our resummed cross sections satisfy the test of stability under 
variation of the hard scale The resummed results show less variation than the NLO cross 
section. On the other hand, this is not true of the resummation of LSvN. As shown in Table 1 
their qq cross section in the DIS scheme has a //-variation of 21%. For comparison, the NLO 
cross section shows a variation of 9% and our resummed cross section a variation of 6%. 

Scheme dependence is an extra source of theoretical uncertainty, but it should produce 
minimal differences for physical cross sections if these are calculated in an unambiguous re- 
summation approach. To check for possible scheme dependent uncertainty, we perform our 
resummation for the dominant qq channel in both schemes. The cross sections presented in 
Table 2 show that our scheme dependence is insignificant, resulting in a difference of about 
4% for the cross section. 



m (GeV) 


(Xqq{DlS; MS) 


/i/m = 0.5 


fi/m = 1 


fi/m = 2 


150 


NLO 


9.42; 9.68 


9.31; 9.53 


8.57 


8.73 






9.76; 10.16 


9.92; 10.42 


9.31 


9.87 


175 


NLO 


4.46; 4.54 


4.39 


4.43 


4.01 


4.02 




pert 


4.63; 4.78 


4.69 


4.87 


4.37 


4.58 


200 


NLO 


2.20; 2.21 


2.15 


2.14 


1.96 


1.93 




pert 


2.29; 2.34 


2.30 


2.37 


2.14 


2.21 



Table 2. Physical cross sections in pb for the qq channel: DIS versus MS scheme. 



4. Discussion and Conclusions 

Our theoretical analysis shows that perturbative resummation without a model for non- 
perturbative behavior is both possible and advantageous. In perturbative resummation, the 
perturbative region of phase space is separated cleanly from the region of non-perturbative 
behavior. The former is the region where large threshold corrections exponentiate but behave 
in a way that is perturbatively stable. The asymptotic character of the QCD perturbative series, 
including large multiplicative color factors, is fiat, and excursions around the optimum number 
of perturbative terms does not create numerical instabilities or intolerable scale-dependence. 
Infrared renormalons are far away from the stability plateau and, even though their presence 
is essential for defining this plateau, they are of no numerical consequence in the perturbative 
regime. Large color factors, which are multiplicative, enhance the IR renormalon effects and 
contribute significantly to limiting the perturbative regime. 

Our resummed cross sections are about 9% above the NLO cross sections computed with 
the same parton distributions. The scale dependence of our cross section is fairly fiat, resulting 
in a 9 — 10% theoretical uncertainty. This variation is smaller than the corresponding depen- 
dence of the NLO cross section, and it is much smaller than the corresponding dependence 
of the resummed cross section of LSvN. There are other perturbative uncertainties, such as 
dependence on parton distributions and factorization scheme, which affect our cross section 
minimally, at level of 4% or less. These variations are strongly correlated, so we opt not to 
add them in estimating the theoretical uncertainty. Commenting briefly on recent papersB, in 
which the authors state that the increase in cross section they find with their resummation 
method is of the order of 1% over NLO, we stand firmly by our resuhs. Asweexplained0there 
is no actual factorial growth in our expressions, of the type they suggest, precisely because we 
stay away from infrared renormalons in the exponent by optimizing the number of perturba- 
tive terms, and because our phase space is constrained within the perturbative regime derived 
by controlling the non-universal subleading logarithms. The numerical difference in the two 
approaches boils down to the treatment of the subleading logarithms, which can easily shift 
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the results by a few percent, if proper care is not taken. Our approach includes the universal 
leading logarithms only while theirs includes non-universal subleading structures which pro- 
duce the suppression they find. In our opinion, their treatment of the subleading structures 
is not correct. We will present an account of these issues in a future publication. 

Our theoretical analysis and the stability of our cross sections under n variation provide 
confidence that our perturbative resummation procedure yields an accurate calculation of the 
inclusive top quark cross section at Tevatron energies and exhausts present understanding 
of the perturbative content of the theory. Our prediction agrees with data, within the large 
experimental uncertainties. 

This work was supported by the US Department of Energy, Division of High Energy 
Physics, Contract W-31-109-ENG-38. 
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